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摘要 : [845] 利用 权重 基因 共 表 达 网 络 分 i ( weighted gene co-expression network analysis, 
WGCNA) IRK & 42:33. Aedes aegypti 不 同 组 织 基因 共 表 达 模 式 。[【 方 法 ] 从 NCBI SRA 数据 库 中 选 
择 埃及 伊 蚊 不 同 组 织 的 转录 组 数据 中 具 代 表 ， dem 9 FFIR (XE IE y, 5c Ak f fo Ris OE NSCAR F A 
JR fe Bp SE , 雄 成 蚊 的 前 足 、 中 足 、 后 足 和 腹部 末端 ) 的 双 端 测序 数据 ;经 过 缺失 值 移 除 以 及 方差 计算 
后 ,筛选 出 方差 最 大 的 5 000 个 基因 ,利用 民 软 件 中 WGCNA 包 建 立 埃及 伊 蚊 成 蚊 不 同 组 织 的 基因 
共 表 达 网 络 并 划分 模块 ;然后 利用 clusterProfiler 包 对 组 织 特 异性 模块 内 的 基因 进行 GO ( Gene 
Ontology ) 和 E ar Encyclopedia of Genes and Genomes) 富 集 分 析 , 并 用 Cytoscape 软件 中 的 
CytoHubba 插件 筛选 共 表 达 模 块 内 的 hub XXIX. [ 4 iesu ed 同 组 织 中 共 鉴 定 出 11 
i I. ERA Rp SE TFAAR SOR LSU RaR PEERI 个 特异 
ERZ EIE, PEPER PARAR, 6 o 隆 表 达 模 块 内 基因 有 功 
能 注释 到 组 织 生物 学 功能 ;其 中 ,此 蚊 触角 特异 性 green 模块 内 基因 具有 气味 结合 和 嗅觉 < 
AE xl] 8E S PESSCUR TE ALTE purple 模块 内 基因 具有 丝氨酸 型 肽 链 内 切 酶 活性 和 丝氨酸 水 解 酶 活性 等 功 
能 ; 雄 蚊 脑 特异 性 blue 模块 内 基因 在 生物 学 过 程 调节 、 信 号 转 导 和 神经 系统 过 程 等 生物 学 过 程 中 
发 挥 主要 作用 。 利 用 CytoHubba 进一步 鉴定 出 所 选 组 织 特 异性 共 表 达 模 块 中 具有 高 连通 性 的 hub 
基因 ,包括 AAEL010426, AAEL002896, AAEL002600, AAEL000961, AAEL007784 和 AAELOO06429 , 
【结论 ] 本 研究 依据 埃及 伊 蚁 不 同 组 织 转录 组 数据 ,利用 WGCNA 方法 发 现 了 许多 重要 的 基因 共 表 
达 模 块 。 本 研究 的 结果 为 豚 吕 基因 共 表 达 模 式 分 析 提 供 新 思路 和 方法 基础 ,对 探究 蚊 忠 不 同 组织 
特有 的 基因 资源 信息 以 及 功能 基因 生物 信息 学 研究 有 参考 价值 。 
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Abstract: [ ^im] To explore the co-expression patterns of genes in different tissues of Aedes aegypti by 
using weighted gene co-expression network analysis ( WGCNA). [Methods] Paired sequencing data of 
representative nine adult tissues ( antenna and brain of female and male, proboscis, maxillary palp and 
ovary of female, and foreleg, midleg, hindleg and abdominal tip of male) were selected from the 
transcriptome data of different tissues of A. aegypti in the NCBI SRA database. After removing missing 


values and calculating variance, 5 000 genes with the highest variance were chosen, and the WGCNA 
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package in R software was used to establish the gene co-expression networks in different adult tissues of 
A. aegypti and divide modules. Then the clusterProfiler package was used to conduct GO ( Gene 
Ontology) and KEGG ( Kyoto Encyclopedia of Genes and Genomes) enrichment analysis of the genes in 
the tissue-specific modules, and CytoHubba plug-in in Cytoscape software was adopted to screen the hub 
genes in co-expression modules. [Results] A total of 11 gene co-expression modules in different adult 
tissues of A. aegypti were identified, and one specific expression module was identified in every tissue of 
the antenna, proboscis, ovary and maxillary palp of female adult, and the brain and abdominal tip of 
male adult. There was no specific expression module in tissues of the foreleg, midleg and hindleg of male 
adults. Genes in six tissue-specific expression modules were annotated to biological functions of the 
corresponding tissues. Among them, genes in the female antenna specific green module have functions 
such as odor binding and olfactory receptor activity. Genes in the female proboscis specific purple module 
have functions such as serine-type endopeptidase activity and serine hydrolase activity. Genes in the male 
brain specific blue module play a major role in regulation of biological process, signal transduction, and 
nervous system processes. The hub genes with high connectivity in selected tissue-specific modules, 
including AAELO10426, AAELO02896, AAELO02600, AAELOO00961, AAELO07784 and AAELO006429 , 
were identified by using CytoHubba. [Conclusion] In this study, based on the transcriptome data of 
different tissues of A. aegypti we discovered many important co-expressed gene modules by using the 
WGCNA method. The results provide new ideas and a method basis for the analysis of mosquito gene co- 


expression patterns and have reference value for exploring the unique gene resource information and 


bioinformatical studies of functional genes in different tissues of mosquitoes. 


Key words: Aedes aegypti; WGCNA; co-expression pattern; transcriptome; hub gene 


截至 2017 年 ,全 球 已 知 蚊虫 约 3 573 种 ,其 中 
我 国 约 有 419 fh TCK, 2018 ) 。 它 们 分 布 
广 , 楷 殖 快 。 绝 大 多 数 种 类 的 肉 蚊 需要 吸食 血液 作 
为 补充 营养 以 使 卵 正常 发 育 。 因 其 吸血 行为 ,蚊虫 
传播 许多 危险 的 人 类 疾病 的 病原 体 ,是 多 种 恶性 传 
染病 的 传播 媒介 (Tolle, 2009) 。 蚊 媒 病 是 人 类 健康 
的 主要 威胁 ,每 年 全 球 有 7 亿 人 受 感染 ,造成 70 多 
万 人 死亡 。 近 年 来 ,登革热 发 病 率 大 幅度 增长 ,全 球 
约 有 一 半 人 口 面 临 登革热 的 威胁 (WHO，2019) 
(https: // www. who. int/zh/ news-room/fact-sheets/ detail/ 





























dengue-and-severe-dengue) , R X Jt. Aedes aegypti 
( Linnaeus, 1762) J& XE E RA Se 0] EREE ES, E] 
时 也 携带 基 孔 肯 雅 热 . 黄 热 病 以 及 寨 卡 病毒 病 等 传 
染病 病原 。 这 些 疾 病 每 年 在 全 世界 造成 数 亿 人 感 
We ,导致 数 万 人 死亡 ,并 呈 逐 年 上 升 趋势 (Akbari et 
al., 2013) 。 媒 介 坚 忠 的 控制 是 蚊 媒 病 控制 的 关键 ， 
已 成 为 当今 社会 所 面临 的 叹 待 解决 的 重大 问题 之 
一 。 蚊 虫 的 不 同 组 织 承 担 不 同 的 生理 功能 ,掌握 它 
们 的 基因 表达 差异 有 利于 认 知 它们 的 生物 学 功能 。 
蚊虫 的 遗传 控制 是 蚊虫 潜在 的 重要 控制 措施 之 一 ， 
































控制 的 靶 标 基因 。 

转录 组 是 特定 细胞 或 组 织 在 特定 时 间或 状态 下 
所 有 转录 产物 的 集合 ( Wang et al., 2010) ,通过 对 转 
录 组 的 研究 可 以 揭示 蚊虫 组 织 分 化 过 程 中 的 基因 表 
达 模 式 和 结构 变异 ,并 发 现 新 基因 等 。 近 年 来 , 随 着 
microarray 和 RNA-Seq 等 高 通 量 技术 的 飞速 发 展 与 
成 本 的 降低 ,产生 了 大 规模 蚊虫 组 织 转录 组 学 数据 
( Martínez-Barnetche et al., 2012; Matthews et al., 
2016; Chen et al., 2017) 。 传 统 的 转录 组 学 研究 主 
要 是 在 两 类 样本 相互 比较 的 基础 上 ,鉴定 差异 表达 
基因 ,进而 对 差异 基因 进行 注释 与 富 集 分 析 ,并 挖掘 
这 些 基因 的 功能 表现 和 参与 的 通路 (pathway) ,从 而 
解析 相关 的 生物 学 问题 或 分 子 机 制 。 但 这 样 的 研究 
策略 主要 局 限于 对 组 间 差 异 表达 的 基因 进行 研究 ， 
容易 忽略 掉 基 因 之 间 的 协同 表达 或 表达 互 作 。 H 
前 ,基于 全 局 的 权重 基因 共 表 达 网 络 分 析 ( weighted 
gene co-expression network analysis, WGCNA ) 正好 为 
3x 4- RESI fe eroe T RHL WGCNA 分 析 基 于 大 
样本 转录 组 数据 ,根据 基因 表达 水 平 的 不 同 将 基因 
聚 类 ,具有 相似 表达 模式 的 基因 会 被 划分 到 同一 个 

































































弄 清 蚊虫 不 同 组织 在 行使 功能 过 程 中 的 共 表 达 基 因 
以 及 起 重要 作用 的 核心 (hub) 基因 有 利于 寻找 遗传 


模块 ,同一 模块 中 的 基因 常 具有 相似 的 生物 学 功能 
( 杨 宇 昕 等 ,2019)。 作 为 一 种 系统 生物 学 分 析 方 
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ik, WGCNA 可 以 用 来 鉴定 高 度 协 同 变化 的 基因 集 、 
挖掘 与 目标 组 织 或 性 状 相关 的 特异 性 模块 ,并 可 进 
一 步 研究 模块 内 基因 的 互 作 关系 ,以 及 筛选 hub 基 
(Langfelder and Horvath , 2008) , 且 在 描述 解析 复 
杂 表 型 下 的 分 子 作用 机 制 和 网 络 关 系 等 方面 具有 独 
RRA, A WGCNA 方法 和 软件 包 开 发 以 来 , 它 已 
广泛 应 用 于 多 个 研究 领域 ,并 取得 了 一 系列 的 研究 
成 果 ( Greenham et al., 2017; Liu et al., 2018; Deng 
et al., 2019) , 

在 蚊虫 转录 组 研究 领域 , Hickner 等 (2015 ) 用 
WGCNA MAWI T REEN Culex pipiens iii A Y 
导 的 潜在 调控 机 制 ,构建 了 包含 5 个 模块 的 基因 共 
表达 网 络 ,通过 对 网 络 模块 中 基因 的 富 集 分 析 , 鉴 定 
出 了 与 滞 育 诱导 相关 的 多 个 基因 集 ;通过 对 埃及 伊 
蚊 环 境 驱 动 变 异 介 导 的 模块 化 基因 网 络 研究 , Kang 
等 (2018) 揭 示 了 先天 人 免疫 反应 的 关键 转录 调控 因 
子 和 关键 基因 。 这 些 研究 表明 ,将 转录 组 学 数据 与 
WGCNA 方法 相 结 合 , 可 为 探索 蚊虫 生物 学 过 程 相 
关 分 子 机 制 提供 强 有 力 的 工具 。 然 而 ,尚未 有 基于 
WGCNA 方法 对 蚊虫 不 同 组 织 基因 共 表 达 模 式 的 
研究 。 

蚊虫 的 不 同 组 织 在 其 传播 疾病 过 程 中 起 着 重要 
的 作用 ,可 作为 防 控 的 靶 标 。 本 研究 利用 WGCNA 
方法 ,结合 RNA-Seq 转录 组 数据 ,构建 了 埃及 伊 蚊 9 
种 组 织 加 权 基 因 共 表达 网 络 ,挖掘 出 与 组 织 特异 性 
相关 的 共 表 达 模 块 ,研究 模块 基因 功能 ,并 筛选 出 在 
组 织 特异 性 模块 内 起 重要 作用 的 hub 基因 ,揭示 基 
间 的 相互 作用 关系 , 旨 在 为 探究 蚊虫 不 同 组 织 特 
有 的 基因 资源 信息 ,进而 解析 基因 功能 与 组 织 功 能 
的 联系 ,筛选 遗传 控制 的 靶 标 基因 提供 新 的 思路 和 


途径 。 
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1 材料 与 方法 


1.1 转录 组 数据 来 源 

埃及 伊 蚊 不 同 组 织 的 转录 组 数据 来 源 于 NCBI 
(National Center for Biotechnology Information ) 的 
SRA ( Sequence Read Archive) 数据 库 ,获取 码 为 
PRJNA236239( Matthews et al., 2016) 。 选 择 其 中 未 
吸血 状态 下 雌 蚊 和 雄 蚊 9 种 组 织 的 双 端 测序 数据 进 
行 分 析 ,包括 肉 蚊 的 触角 、 脑 ` 喉 .下 频 须 、 卵 梨 , 以 及 
雄 蚊 的 触角 、 脑 ` 前 足 . 中 足 \ 后 足 和 腹部 末端 (被 定 
XN Hilo EL E PH RE ALBUS 3 个 末端 腹 节 ) 等 。 不 
同 组 织 数 据 包 含 的 生物 学 重复 次 数 为 : 雌 蚊 咏 5 次 











生物 学 重复 ,触角 、 下 频 须 4 次 重复 ,卵巢 3 次 重复 ， 
脑 2 次 重复 , 雄 蚊 各 组 织 均 3 次 重复 ,总 计 36 个 样 
本 。WCCNA 分 析 根 据 基因 表达 水 平 的 不 同 将 基因 
聚 类 。 在 转录 组 研究 中 ,可 通过 比 对 到 基因 组 区 域 
或 基因 外 显 子 区 的 测序 序列 (reads) 计数 ( read 
counts ) 来 估计 基因 的 表达 水 平 。 但 reads 计数 除了 
与 基因 的 真实 表达 水 平成 正比 外 ,还 与 基因 的 长 度 
和 测序 深度 成 正 相 关 , 因 此 使 得 不 同 基 因 的 表达 水 
平 不 具有 可 比 性 。 为 了 解决 这 个 问题 ,人 们 引入 了 
FPKM( fragments per kilobase million ) , RPKM ( reads 
per kilobase million ) 和 TPM ( transcripts per million) 
的 概念 。 其 中 TPM 是 每 百 万 条 reads 的 转录 本 ,其 
同时 考虑 了 测序 深度 和 基因 长 度 对 reads 计数 的 影 
响 ,是 目前 常用 的 衡量 基因 表达 水 平 的 方法 之 一 。 
本 研究 中 埃及 伊 蚊 不 同 组 织 基 因 的 TPM 值 采用 
Matthews 等 (2016 ) 文献 中 Additional File 6 中 的 


+ HH 
结果 。 


1.2 基因 共 表 达 网 络 构建 

采用 TPM 值 作 为 标准 化 的 基因 表达 水 平定 量 ， 
本 研究 共 纳 入 36 个 样本 。 在 应 用 WGCNA 前 ,基于 
公式 log, (TPM + 1) 对 数据 进行 转换 ( 赵 志 洪 ， 
2017) 。 经 过 缺失 值 移 除 以 及 方差 计算 后 ,筛选 出 
方差 最 大 的 5 000 个 基因 ,使 用 R(v.3.6.1) 软 件 中 
(https: // cran. r-project. org/ bin/ windows/ base/old/ ) 
的 WGCNA(v.1.67) 包 来 构建 加 权 基 因 共 表达 网 络 
并 划分 模块 。 首先 利用 网 数 pickSoftThreshold 
( Langfelder and Horvath, 2008 ) 计算 权重 值 , 通 过 选 
择 合适 的 权重 值 ,使 网 络 符合 无 尺度 网 络 分 布 。 

根据 上 一 步 分 析 结 果 , 选 择 最 优 的 软 阐 值 。 使 
H A x ed 28 F4 S£ ER žit blockwiseModules ( Langfelder 
and Horvath, 2008) 构建 共 表 达 网 络 , 除 参 数 
minModuleSize =50 外 ,其 余 参 数 按照 默认 设置 。 
1.3 基因 表达 组 织 特异 性 模块 鉴定 

一 组 表达 模式 高 度 相似 的 基因 常 聚集 在 同一 模 
块 ,而 具有 相似 表达 模式 的 基因 被 认为 具有 类 似 的 
生物 学 功能 ( 杨 宇 昕 等 , 2019 ) 。 一 般 而 言 , 如 果 某 
模块 与 样品 的 某 种 特定 表 型 的 相关 性 显著 高 于 其 他 
模块 ,说 明 这 一 个 模块 可 能 与 该 表 型 存在 强 的 关联 ， 
为 特异 性 模块 。 本 研究 由 于 只 有 埃及 伊 蚊 不 同 组 织 
的 转录 组 数据 ,并 没有 其 他 表 型 数据 ,在 此 直接 将 不 
同 组 织 作 为 研究 的 表 型 性 状 ,以 分 析 模 块 与 组 织 的 
联系 。 首 先 建立 一 个 行为 36 个 组 织 样本 , 列 为 11 
个 不 同 取样 组 织 名 称 的 抢 阵 ,者 行 与 列 为 同一 组 织 ， 
BIS 1 ,不 是 同一 组 织 即 为 0。 计 算 每 个 模块 特征 向 
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量 基 因 (module eigengene, ME) 和 这 个 和 矩阵 之 间 的 
相关 系数 ,该 数值 越 接 近 1 ,表示 模块 与 样品 表 型 的 
正 相 关 性 越 强 ; 越 接近 - 1 ,表示 模块 与 样品 表 型 负 
相关 性 越 强 。 本 研究 中 模块 与 组 织 高 度 正 相关 的 筛 
选 标准 为 :模块 与 组 织 相 关系 数 大 于 0. 65 , 且 统 计 
显著 性 己 <0.05。 选 择 与 组 织 高 度 正 相关 的 模块 ， 
作为 组 织 特异 性 模块 进行 后 续 分 析 ( Downs et al., 
2013 ) 。 
1.4 组织 特异 性 模块 基因 的 GO KEGG 富 集 
分 析 

首先 ,使 用 Bioconductor 中 的 AnnotationHub ( v. 
2.18.0) 包 获取 埃及 伊 蚊 的 注释 信息 ,下 载 埃及 伊 
Iit OrgDb 数据 库 ( AH76482, v.2.1) 中 的 注释 数据 ， 
再 利用 Bioconductor 中 的 clusterProfiler ( v. 3. 14. 3) 
包 对 模块 内 基因 进行 CO( Gene Ontology ) 和 KEGG 
( Kyoto Encyclopedia of Genes and Genomes ) 富 集 分 
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析 (Yu et al., 2012), 一般 情况 下 ,P <0.01,g < 
0. 05 , 则 为 宣 集 功能 项 或 通路 。 
1.5 模块 基因 互 作 网 络 构建 筑 选 hub 基因 

导出 WGCNA 共 表 达 网 络 分 析 结 果 到 
Cytoscape(v.3.7.1) 软 件 (Su et al., 2014) 中 ,对 基 
因 调控 网 络 进 行 可 视 化 。 通 过 CytoHubba 插件 的 12 
种 算法 ,得 到 连接 度 在 前 10 的 基因 ,选择 在 不 同 算 
法 中 重复 出 现 次 数 大 于 5 的 基因 ,在 本 研究 中 被 定 
义 为 hub 基因 。 






































2 结果 

2.1 权重 基因 共 表 达 网 络 构建 

根据 基因 表达 水 平 对 样本 进行 聚 类 。 如 图 1 所 
基于 平均 (average) 距 离 和 层次 聚 类 算法 聚 类 结 
果 表 明 ,基因 表达 模式 除了 在 雄 成 蚊 前 足 和 中 足 中 
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基于 基因 表达 水 平 的 埃及 伊 蚊 成 蚊 不 同 组 织 聚 类 
Clustering of different adult tissues of Aedes aegypti based on gene expression levels 
Fe Ov; 雌 蚊 卵 梨 Female ovary; Ma, Br; 雄 蚊 脑 Male brain; Fe, Br: WEZ Female brain; Ma, An: 雄 蚊 触角 Ma 
Female antenna; Ma, At; 雄 蚊 腹部 末端 Male abdominal tip; Ma HL; 雄 蚊 后 足 Male hindleg; Ma FL; 雄 蚊 前 足 Male foreleg; Ma ML; HEP E 
Male midleg; Fe Os; WE Female proboscis; Fe Pa; 肉 蚊 下 显 须 Female maxillary palp。 根 据 基 
样本 进行 聚 类 。 横 轴 表 示 埃及 伊 蚊 成 蚊 不 同 组 织 ,代号 未 位 的 数字 表示 该 组 织 的 重复 样本 ; 纵 轴 代 表 基 
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,采用 平均 距离 和 层次 聚 类 算法 对 
间 的 聚 类 高 度 。 


























According to gene 





expression level, average distance and hierarchical clustering algorithm were used to cluster the samples. The abscissa represents the different adult tissues 


of A. aegypti, the last digit of the code represents a duplicate sample of the tissue, and the ordinate represents the cluster height of genes. 
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出 现 了 一 定 的 散乱 之 外 ,测序 所 用 的 埃及 伊 蚊 成 蚊 
组 织 大 多 具有 很 高 的 重复 性 。 雌 / 雄 成 蚊 的 触角 A 
中 基因 表达 模式 相近 , 聚 类 在 一 起 ; 雄 成 蚊 前 足 、 中 
足 及 后 足 中 基因 表达 模式 相近 , 聚 类 在 一 起 ; 雌 成 蚊 
喉 与 下 咯 须 中 基因 表达 模式 相近 , 聚 类 在 一 起 。 对 
AX BUE CB) 进行 筛选 ,最 终 选 取 B 值 为 12。 通 过 
WGCNA 方法 成 功 构建 埃及 伊 蚊 成 蚊 不 同 组 织 加 权 
基因 共 表 达 网 络 ,并 将 网 络 划 分 为 11 个 模块 (图 2: 
A)。 图 2(A) 上 半 部 分 的 聚 类 树 中 每 个 校 代表 一 个 
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模块 颜色 $ 
Module colors 


Gene number 


模块 ,每 个 枝 上 的 叶 代表 一 个 基因 ,图 2(A) 下 半 部 
分 的 模块 颜色 中 每 个 颜色 代表 一 个 模块 。 共 表达 模 
块 中 基因 的 具体 数量 如 图 2(B) 所 示 。 模 块 中 基因 
个 数 123 ~1 423 不 等 , 共 含 有 5 000 个 基因 。 其 中 
turquoise 模块 包含 的 基因 数目 最 多 , 达 1 423; 
greenyellow 模块 内 的 基因 数量 最 少 ,为 123 ; grey K 
示 未 划分 到 模块 中 的 基因 ,图 2(A) 中 grey 模块 基 
本 不 可 见 , 因 未 划分 到 模块 的 基因 仅 有 7 个 。 
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图 2 埃及 伊 蚊 成 蚊 不 同 组 织 基因 共 表 达 模 块 
Fig. 2 Gene co-expression modules in different adult tissues of Aedes aegypti 
A: 基因 聚 类 数 和 模块 分 制 。 上 半 部 分 表示 基因 聚 类 所 得 聚 类 树 , 纵 坐 标 代表 各 基因 间 的 聚 类 距离 ;下 半 部 分 表示 按 树 的 分 梳 切 割 得 到 的 模 
块 ,相同 模 块 用 同一 颜色 表示 。 Clustering dendrogram of genes and module splitting. The upper part represents dendrogram tree obtained by gene 




















hierarchical clustering, and the ordinate represents the clustering distance between genes. The lower half represents the module cut by the branches of the 
tree, and the same module is represented by the same color. B: 共 表 达 模 块 中 基因 数量 分 布 。 横 坐标 代表 模块 , 纵 坐 标 代表 模块 中 的 基因 数目 ,不 
同 的 颜色 与 模块 名 称 对 应 ,每 个 柱子 上 的 数字 表示 该 模块 内 的 基因 数目 。Distribution of gene number. The abscissa represents the module, and the 


ordinate represents the number of genes in the module. Different colors correspond to the module name. The number on each column represents the 





number of genes in the module. 


2.2 ”组织 特异 性 模块 分 析 结 果 

通过 往 选 发 现在 11 个 模块 中 有 6 个 为 高 度 的 
组 织 特异 性 模块 (r >0.65, P «0.05) (R13). fle 
成 蚊 触 角 咏 、 卵 和 偶 、 下 颗 须 以 及 雄 成 蚊 脑 和 腹部 末 
端 组 织 中 各 鉴定 出 1 个 特异 性 模块 。green 模块 与 
雌 成 蚊 触 角 存在 高 度 相 关 性 (r =0. 69, P =4E -06), 
purple 模块 与 肉 成 蚊 唆 存在 高 度 相关 性 (r= 0. 8， 
P=4E -09) , turquoise 模块 与 肉 成 蚊 卵 集 存在 高 度 
正 相 关 性 (r =0.96, P=4E - 20) ,red 模块 与 肉 成 蚊 
下 显 须 存在 高 度 正 相关 性 (r=0.69, P -4E -06)。 
与 雄 成 蚊 脑 高 度 正 相 关 的 模块 为 blue 模块 (+ = 
0.7, P=2E -06) ,与 雄 成 蚊 腹部 末端 高 度 正 相关 








足 . 中 足 \ 后 足 存在 较 高 相关 性 。 推 测 可 能 是 这 3 个 
组 织 中 的 基因 主要 进行 基本 的 生命 活动 ,因此 没有 
特异 性 的 模块 与 之 对 应 。 
2.3 组 织 特异 性 模块 GO 和 KEGG 富 集 分 析 

组 织 特 异性 模块 可 能 参与 组 织 特 有 的 生物 学 过 
FE ,为 进一步 解析 组 织 特异 性 模块 内 基因 的 生物 学 
功能 ,并 与 相应 组 织 功能 联系 ,同时 也 验证 共 表达 网 
络 构建 和 组 织 特异 性 模块 的 可 靠 性 ,对 这 6 个 特异 
性 模块 进行 GO 和 KEGG ARIT GO 富 集 分 析 
结果 表明 ( 表 1) ,与 肉 成 蚊 触角 高 度 正 相关 的 green 
模块 在 纤毛 组 装 、 基 于 微 管 运动 和 细胞 对 刺激 的 反 














的 模块 为 magenta 模块 (r =0.96, P =2E -37), 6 
个 组 织 特异 性 模块 中 没有 任何 一 个 模块 和 雄 成 蚊 前 


应 等 生物 学 过 程 中 发 挥 主要 作用 ,其 主要 的 细胞 组 
分 为 微 管 相关 复合 物 、 细 胞 外 围 . 质 膜 等 ,具有 气味 
结合 嗅觉 受 体 活性 . 览 膜 信号 受 体 活 性 等 生物 学 功 
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图 3 埃及 伊 蚊 成 蚊 不 同 组织 基 因 共 表达 模块 与 性 状 关联 热 图 





Fig. 3 Heat map of module-trait relationship of different adult tissues in Aedes aegypti 
Fe An; 肉 蚊 触角 Female antenna; Fe Br; MERCAR Female brain; Fe Os; WEE Female proboscis; Fe Ov: 肉 蚊 卵 梨 Female ovary; Fe Pa; WHfEd 
“ji Female maxillary palp; Ma, An; 雄 蚊 触 角 Male antenna; Ma. At; 雄 蚊 腹部 末端 Male abdominal tip; Ma. Br; 雄 蚊 脑 Male brain; Ma, FL; 雄 蚊 
前 足 Male foreleg; Ma ML: 雄 蚊 中 足 Male midleg; Ma HL; 雄 蚊 后 足 Male hindleg. MEgreenyellow: 黄 绿 色 模 块 特征 向 量 基因 Module eigengene 




















of green yellow module; MEturquoise: 蓝 绿色 模块 特征 向 量 基 因 Module eigengene of turquoise module; MEmagenta: 品 红 色 模 块 特 征 向 量 基 因 
Module eigengene of magenta module; MEpurple; 紫色 模块 特征 向 量 基 因 Module eigengene of purple module; MEred: 红色 模块 特征 向 量 基 因 


























Module eigengene of red module; MEblack : 色 模 块 特征 向 量 基因 Module eigengene of black module; MEbrown: 棕色 模块 特征 向 量 基因 Module 
eigengene of brown module; MEyellow: 黄色 模块 特征 向 量 基因 Module eigengene of yellow module; MEblue: 蓝 色 模块 特征 向 量 基因 Module 
eigengene of blue module; MEgreen: 绿色 模块 特征 向 量 基因 Module eigengene of green module; MEpink: 粉红 色 模 块 特征 向 量 基 因 Module 
eigengene of pink module; MEgrey: 灰色 模块 特征 向 量 基 因 Module eigengene of grey module. 该 图 基于 R 语言 所 得 , 横 轴 表示 埃及 伊 蚊 不 同 组 织 ， 
纵 轴 表示 每 一 个 模块 的 特征 向 量 。 括 号 内 数值 表示 模块 与 组 织 相 关系 数 和 统计 显著 性 。 相 关系 数 越 接近 1 ,表示 模块 与 样品 正 相 关 性 越 强 ; 
越 接 近 - 1 ,表示 模块 与 样品 负 相 关 性 越 强 。Based on the R language, the horizontal axis represents the different tissues of A. aegypti, and the vertical 


axis represents the module eigengene of each module. The values in parentheses represent the correlation coefficient and statistical significance between the 



















































































module and the tissues. The closer the correlation coefficient is to 1, the stronger the positive correlation between the module and the sample; the closer 


to —1, the stronger the negative correlation between the module and the sample. 


能 ;与 雄 成 蚊 脑 高 度 正 相 关 的 blue 模块 在 生物 学 过 
程 调节 、 信 和 号 转 导 、 神 经 系统 过 程 等 生物 学 过 程 中 发 
挥 主 要 作用 ,其 主要 的 细胞 组 分 为 质 膜 、 突 触 、 突 触 
后 膜 等 ,具有 转运 活性 、 神 经 递 质 受 体 活性 、 铁 离子 
结合 等 功能 ;与 肉 成 蚊 唆 高 度 正 相 关 的 purple 模块 
具有 丝氨酸 型 肽 链 内 切 栈 活性、 丝氨酸 水 解 酶 活性 
等 功能 ;与 雌 成 蚊 卵 梨 高 度 正 相关 的 turquoise 模块 
Tr DNA 复制 、 细 胞 分 裂 、 细 胞 周期 等 生物 学 过 程 中 
发 挥 主要 作用 ,其 主要 的 细胞 组 分 为 核 小 体 .染色 体 
等 ,具有 DNA 结合 、 锌 离子 结合 、 微 管 结合 、 解 旋 酶 
活性 等 功能 ; 与 肉 蚊 下 频 须 高 度 正 相关 的 red 模块 
在 有 机 酸 代 谢 、 羚 酸 代 谢 、 细 胞 氨基 酸 代 谢 等 生物 学 
过 程 中 发 挥 作用 ,其 细胞 组 分 为 细胞 外 区 域 , 具 有 丝 
氮 酸 型 肽 链 内 切 酶 活性 碳水 化 合 物 结合 等 功能 ;与 
雄 成 蚊 腹 部 末端 高 度 正 相关 的 magenta Eit fr 2 We 
代谢 .草酸 代谢 ` 有 机 酸 代谢 等 生物 学 过 程 中 发 挥 主 
要 作用 。 

KEGG 富 集 分 析 结 果 表 明 ( 表 2) ,与 肉 成 蚊 触 














角 高 度 正 相关 的 green 模块 显著 富 集 在 药物 代谢 - 细 
胞 色素 P450 ,不 饱和 脂肪 酸 的 生物 合成 以 及 谷 胱 甘 
肽 代谢 等 通路 ;与 雄 成 蚊 脑 高 度 正 相关 的 blue 模块 
显著 富 集 在 刺激 神经 组 织 中 的 交互 、 光 传导 以 及 其 
他 多 糖 降解 通路 ; ts ME po coc BD E. eg BE IE IH OR I 
turquoise 模块 显著 富 集 在 DNA 复制 、 错 配 修 复 、 戊 
糖 和 和 葡萄糖 醛 酸 盐 的 相互 转化 等 通路 ;与 雌 成 蚊 下 
颗 须 高 度 正 相关 的 red. 模块 显著 富 集 在 甘氨酸 、 丝 
氮 酸 和 苏 氮 酸 的 代谢 通路 以 及 精 所 酸 和 有 情 氨 酸 代谢 
通路 ;magenta 和 purple 模块 未 发 现 显著 富 集 通路 。 
2.4 组 织 特异 性 模块 hub 基因 识别 

整合 CytoHubba 插件 12 种 不 同 算法 的 结果 , 取 
重复 出 现 次 数 大 于 5 次 的 基因 , 共 得 到 47 个 hub A 
(green, blue, purple, turquoise, red 和 magenta 模 
块 分 别 有 9, 5, 6, 9, 9 $09 个 ) (XX 3). green, 
blue, purple, turquoise, red 和 magenta 模块 中 连通 
性 最 高 的 基因 分 别 为 AAELO010426, AAELOO2896 , 
AAELO02600 , 4457000961, AAELOO778A 和 AAELOO6429 , 
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表 1 埃及 伊 蚁 成 蚁 组 织 特 异性 基因 共 表达 模块 的 GO 富 集 分 析 
Table 1 GO enrichment analysis of tissue specific gene co-expression module in Aedes aegypti adults 

本 体 描述 q f& 模块 特异 性 相关 组 织 

Ontology Description q-value Module Specifically related tissue 

BP 纤毛 组 装 Cilium assembly 1.35E -03 Green Efi ffi Female antenna 

BP 基于 微 管 运动 Microtubule-based movement 2.47E -03 Green Efi ffi Female antenna 

BP 细胞 对 刺激 的 反应 Cellular response to stimulus 4.23E - 17 Green Efi ff] Female antenna 

CC 质 腊 Plasma membrane 2.25E - 18 Green 蚊 触 角 Female antenna 

CC 细胞 外 围 Cell periphery 5.56E - 18 Green 上 蚊 触 角 Female antenna 

CC 微 管 相 关 复 合 物 Microtubule associated complex 1.38E -03 Green Efi ff Female antenna 

MF 气味 结合 Odorant binding 1.67E -45 Green 蚊 触 角 Female antenna 

MF 嗅觉 受 体 活性 Olfactory receptor activity 7.97E -43 Green Erf ffi Female antenna 

MF 跨 膜 信号 受 体 活性 Transmembrane signaling receptor activity 5.78E -31 Green 蚊 触 角 Female antenna 

BP 生物 学 过 程 调 节 Regulation of biological process 2.44E -05 Blue 雄 收 脑 Male brain 

BP 信和 号 转 导 Signal transduction 7.80E -08 Blue 雄 蚊 脑 Male brain 

BP 神经 系统 过 程 Nervous system process 6.58E -03 Blue 雄 蚊 脑 Male brain 

CC 质 腊 Plasma membrane 1.78E -03 Blue 雄 蚊 脑 Male brain 

CC 突 触 Synapse 1.25E -04 Blue HERUR Male brain 

CC 突 触 后 膜 Postsynaptic membrane 1.25E -04 Blue 雄 蚊 脑 Male brain 

MF 转运 活性 Transporter activity 8.97E - 06 Blue 雄 蚊 脑 Male brain 

MF 神经 递 质 受 体 活性 Neurotransmitter receptor activity 6.60E -04 Blue 雄 蚊 脑 Male brain 

MF 铁 离 子 结 合 Iron ion binding 3.73E - 03 Blue 雄 蚊 脑 Male brain 

MF 丝氨酸 型 肽 链 内 切 酶 活性 Serine-type endopeptidase activity 4.78E -10 Purple [Eti Female proboscis 

MF 丝氨酸 水 解 酶 活性 Serine hydrolase activity 4.78E - 10 Purple [ESI Female proboscis 

BP DNA 复制 DNA replication 2.53E - 05 Turquoise 惟 蚊 卵 策 Female ovary 

BP 细胞 分 裂 Cell division 8.52E -04 Turquoise Ec Bg. Female ovary 

BP 细胞 周期 Cell cycle 1.48E -03 Turquoise [Edi BE S& Female ovary 

CC 核 小 体 Nucleosome 8.46E -03 Turquoise HEB gH Female ovary 

CC 染色 体 Chromosome 5.14E -06 Turquoise Edi BE S& Female ovary 

MF DNA 结合 DNA binding 2.02E -04 Turquoise Esc p. Female ovary 

MF 锌 离子 结合 Zinc ion binding 5.15E - 03 Turquoise KE S& Female ovary 

MF 微 管 结合 Microtubule binding 4.33E -03 Turquoise Em B. Female ovary 

MF 解 旋 酶 活性 Helicase activity 4.33E -03 Turquoise Eb gH Female ovary 

BP 有 机 酸 代 谢 过 程 Organic acid metabolic process 2.02E -03 Red ME FS Female maxillary palp 

BP 细胞 氨基 酸 代 谢 过 程 Cellular amino acid metabolic process 5.61E -03 Red MEt P SS Female maxillary palp 

BP 羧 酸 代谢 过 程 Carboxylic acid metabolic process 5.61E -03 Red WEB FS Female maxillary palp 

CC 细胞 外 区 域 Extracellular region part 3.74 - 04 Red Jio F SS Female maxillary palp 

MF 丝氨酸 型 肽 链 内 切 酶 活性 Serine-type endopeptidase activity 1.45E - 10 Red MEt P SS Female maxillary palp 

MF 碳水 化 合 物 结合 Carbohydrate binding 5.55E -03 Red MEOT S Zi Female maxillary palp 

BP 羧 酸 代谢 过 程 Carboxylic acid metabolic process 3.45E -03 Magenta 雄 蚊 腹部 末端 Male abdominal tip 

BP 有 机 酸 代 谢 过 程 Organic acid metabolic process 3.45E - 03 Magenta 雄 蚊 腹部 末端 Male abdominal tip 

BP 草酸 代谢 过 程 Oxoacid metabolic process 3.45E -03 Magenta 雄 蚊 腹部 末端 Male abdominal tip 
BP: 生物 学 过 程 Biological process; CC; 细胞 组 分 Cellular component; MF: 分 子 功能 Molecular function. q {É q-value: 校正 后 的 P 值 Corrected P- 
value 


基于 WGCNA 共 表 达 网 络 分 析 结 果 , 对 组 织 特 





和 4451013563; blue 模块 的 基因 在 肉 / 雄 成 蚊 脑 中 





异性 模块 内 连通 性 最 高 的 核心 基因 进行 可 视 化 展 
示 , 因 与 核心 基因 互 作 的 基因 较 多 , 仅 对 权重 值 前 
50 个 有 关联 的 基因 进行 可 视 化 (图 4)。 结 果 发 现 ， 
green 模块 的 基因 在 雌 / 雄 成 蚊 触 角 中 高 表达 , 与 核 
心 基因 AAELO10426 具有 较 高 互 作 网 络 关系 的 基因 为 
AAELO17505, AAELO17129, AAELO17104, AAELO07147 





高 表达 ,与 核心 基因 AAELO02896 具有 较 高 互 作 网 
络 关系 的 基因 为 44EL012541，AAEL004246， 
AAEL007559, AAEL003472 和 AAELOO2220; purple 
模块 的 基因 在 肉 成 蚊 唆 中 表达 最 高 ,与 核心 基因 
AAEL002600 具有 和 较 高 互 作 网 络 关 系 的 基因 为 
AAELO14382 , AAELO17265 , AAELO11610, AAELO13032 
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R2 埃及 伊 蚊 成 蚊 组 织 特异 性 基因 共 表 达 模 块 的 KEGG 富 集 分 析 


Table 2 KEGG enrichment analysis of tissue specific gene co-expression module in Aedes aegypti adults 


描述 


Description 





药物 代谢 -细胞 色素 P450 Drug metabolism-cytochrome P450 
不 饱和 脂肪 酸 的 生物 合成 Biosynthesis of unsaturated fatty acids 
A EH KIA Glutathione metabolism 

光 传 导 Phototransduction 

刺激 神经 组 织 中 的 交互 Neuroactive ligand-receptor interaction 
糖 降 解 Other glycan degradation 

DNA 复制 DNA replication 

背 腹 轴 形 成 Dorso-ventral axis formation 


错 配 修 复 Mismatch repair 














= 

















T CHAUSSURE VI Arginine and proline metabolism 
q fü q-value; 校正 后 的 P 值 Corrected P-value. 


戊 糖 和 葡萄 糖 醋 酸 盐 的 相互 转化 Pentose and glucuronate interconversions 
甘氨酸 丝氨酸 和 苏 氨 酸 的 代谢 Glycine, serine and threonine metabolism 






































4 值 模块 特异 性 相关 组 织 

q-value Module Specifically related tissue 
2.05E -03 Green 雌 蚊 触角 Female antenna 
2.05E -03 Green 雌 蚊 触角 Female antenna 
4. 56E -03 Green 雌 蚊 触角 Female antenna 
6.32E -06 Blue 雄 蚊 脑 Male brain 
3.39E -05 Blue 雄 蚊 脑 Male brain 
3. 18E -04 Blue 雄 蚊 脑 Male brain 
5.04E -07 Turquoise [Edi BE S& Female ovary 
3.93E -03 Turquoise 惟 蚊 卵 策 Female ovary 
5.49E -03 Turquoise 准 蚊 卵 梨 Female ovary 
7.13E -03 Turquoise 准 蚊 卵 策 Female ovary 
5.02E -04 Red ME F SS Female maxillary palp 
9. 49E -04 Red MEt F SS Female maxillary palp 








RI 埃及 伊 蚊 成 蚊 组 织 特异 性 基因 共 表 达 模 块 内 的 hub 基因 


Table3 Hub genes within a tissue-specific module in Aedes aegypti adults 


模块 Module 2H 2H Tissue 


核心 基因 Hub gene 
































Jis f AAELO010426 , 
Green 

Female antenna 
Blue 雄 蚊 脑 Male brain 
Purple WEE Female proboscis 

MEIK AAELO006726 , 
"Turquoise 

Female ovary 

MFA 
Red 

Female maxillary palp 

MEIBOR EBA Yi AAELO06427 , 
Magenta 








Male abdominal tip 


fil AAELO14390 ; turquoise $R ER HiAk DS] TE ME IX IN K 
中 高 表达 ,与 核心 基因 AAEL000961 具有 较 高 互 作 
网 络 关 系 的 基因 为 A4EL000898，AAELO013338， 
AAEL000923 , AAELOO7879 和 4A4EL011516; red 模块 
的 基因 在 上 肉 成 蚊 下 颇 须 、 雄 成 蚊 腹部 末端 较 高 表达 ， 
与 核心 基因 AAELO07784 具有 较 高 互 作 网 络 关 系 的 
基因 为 AAELOO0044, AAEI003114,  AAELOO7083, 
AAELO02624 FI AAELOI2726 ; magenta 模块 的 基因 在 雄 
成 蚊 腹 部 末端 高 表达 ,与 核心 基因 AAELO06429 H 
有 较 高 互 作 网 络 关 系 的 基因 为 AAELOT0040 , 
AAELOO6414 , AAELOO2000 , AAELOO1550 和 AAELOO5107 , 
已 知 AAELOO7373 和 AAELOOS046 可 能 参与 埃 
及 伊 蚊 寻找 宿主 过 程 中 氮 的 检测 ( Matthews. et al., 
2016) 。 基 于 WGCNA 分 析 结 果 , 发 现 AAELOOT373 
被 划分 到 green 模块 ,该 模块 基因 在 肉 / 雄 成 蚊 触 角 
中 高 表达 ;44EL008046 被 划分 到 blue 模块 ,该 模块 








AAEL001224, AAEL005810, AAEL007669, AAELO010214, AAELO12424, 
AAELO13197 , AAELOO01510 , AAELO14089 

AAELO10918, AAELO02896 , AAELOO1696 , AAELOO2713 , AAELO12850 

AAELO00011, AAELO00039 , AAELOO2600 , AAELO05093 , AAELO13033, AAELO13032 
AAELO008595, AAELO09263, AAEL000961, AAELO003226, AAELO004831 , 
AAELO009654 , AAELO15127 , AAELO14276 

AAELO00044, AAEL007784, AAELO11624, AAELOO7033, AAEL000311, AAELO003405, 
AAELO006920 , AAELOO07102 , AAELO13417 

AAELO17027, AAELO006429, AAELO12488, AAELO17567, AAEL004679 , 
AAELO004944 , AAELO06418 , AAELO10040 


基因 在 上 肉 / 雄 成 蚊 脑 中 高 表达 。 以 这 两 个 基因 为 核 
心 构建 局 部 网 络 ( 图 5) ,发 现 green 模块 内 的 基 
AAELO07373 与 离子 通道 型 受 体 基因 AAELOO0047 
气味 受 体 基 因 AAELO15147 和 AAELOOSO99 等 具有 
较 高 互 作 网 络 关系 ; green 模块 内 的 核心 基 
AAELO14089, AAELOO7669, AAELO12424, AAELO10426, 
AAEL001510 和 AAELO13197 也 均 与 其 存在 互 作 。 
blue 模块 内 的 基因 AAEL008046 与 AAEL004246， 
AAEL006825 和 AAELO00504. 等 存在 关联 ;blue 模块 
内 的 核心 基因 AAELOO2896, AAELO10918 和 
AAELO01696 也 与 其 存在 互 作 网 络 关 系 。 














3 讨论 与 结论 


基因 的 表达 存在 相关 性 ,功能 相关 的 一 组 基因 
第 协调 一 致 ,共同 表达 。 对 蚊虫 不 同 组 织 的 基因 共 
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刘 RAE. 埃及 伊 蚊 不 同 组 织 的 基因 共 表 达 模 式 分 析 1179 
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图 4 埃及 伊 蚊 成 蚊 组 织 特异 性 基因 共 表 达 模 块 内 连通 性 最 高 的 核心 基因 的 基因 网 络 























Fig. 4 Gene networks of the most interconnected hub gene in the tissue-specific gene co-expression module of Aedes aegypti adults 














图 中 绿色 、 蓝 色 、 紫 色 、 蓝 绿色 红色 和 品 红色 线条 分 别 表示 green, blue, purple, turquoise, red 和 magenta 模块 。 每 个 网 络 图 中 位 于 中 心 的 为 该 












































组 织 特异 性 基 共 表达 模块 内 连通 性 最 高 的 核心 基因 ,其 余 带 有 颜色 标记 的 基 表示 与 模块 内 连通 性 最 高 的 核心 基因 具有 和 较 高 互 作 网 络 关系 
的 前 5 个 基因 ,颜色 越 深 表 示 互 作 关 系 越 强 。 网 络 中 只 列 出 与 核心 基因 互 作 权重 值 较 高 的 基因 。The green, blue, purple, turquoise, red, and 


magenta lines in the figure represent the green, blue, purple, turquoise, red, and magenta modules, respectively. At the center of each network diagram 



































is the hub gene with the highest connectivity within the tissue-specific gene co-expression module. The remaining genes with color marks represent the top 


5 genes with a higher interactive network relationship with the hub genes with the highest connectivity in the module. The darker the color, the stronger the 


interactive relationship. In the network, only the genes with high weights interacted with the hub genes were listed. 
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图 5 埃及 伊 蚊 AAELOO7373 和 AAELOOS046 的 基因 网 络 
Fig. 5 Gene networks of AAELO07373 and AAELO008046 of Aedes aegypti 
中 绿色 和 蓝 色 分 别 表示 green 和 blue 模块 。 网 络 中 只 列 出 与 核心 基因 互 作 权重 值 较 高 的 基因 。Green and blue colors in the figure represent 








green and blue modules, respectively. In the network, only the genes with high weights interacted with the hub genes were listed. 





表达 模式 进行 研究 ,可 获得 组 织 特有 的 基因 资源 信 ”中 基因 间 的 相互 作用 模式 ,进而 解析 基因 功能 与 组 
息 ,得 到 组 织 在 行使 功能 过 程 中 的 共 表达 基因 以 及 织 生 物 学 功能 的 联系 ,筛选 遗传 控制 的 训 标 基因 。 
起 重要 作用 的 核心 基因 ,帮助 研究 者 深入 理解 组 织 。 在 大 样本 转录 组 数据 中 ,传统 的 差异 分 析 或 趋势 分 
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析 , 仅 罗列 了 单一 的 基因 列表 ,无 法 得 到 基因 间 的 相 
关 性 ,并 对 基因 进行 有 效 的 分 类 , 亦 不 能 推测 新 的 基 
因 相互 作用 关系 。WGCCNA 分 析 从 系统 角度 出 发 ， 
将 协同 变化 的 基因 进行 聚 类 , 以 网 络 模块 的 形式 呈 
现 基因 的 共 表 达 特 性 。 其 优势 在 于 ,将 海量 的 数据 
进行 归纳 和 整理 ,高 效 研 究 基 因 整 体 表达 规律 ,并 系 
统 地 反馈 样本 中 基因 间 的 相互 作用 模式 ,同时 可 特 
异性 地 筛选 出 与 样本 性 状 相关 的 基因 模块 和 hub 基 
因 。 作 为 一 种 高 效 的 数据 挖掘 手段 ,其 在 多 样本 
RNA-Seq 数据 的 应 用 中 越 来 越 流 行 (Zhao et al., 
2010; 巨 飞 燕 等 , 2019)。WGCNA 分 析 可 为 蚊虫 不 
同 组 织 的 基因 共 表 达 模 式 研究 提供 新 思路 。 

本 研究 利用 WGCNA 方法 对 埃及 伊 蚊 成 蚊 不 同 
组 织 的 基因 共 表 达 模 式 进行 研究 ,构建 了 埃及 伊 蚊 
9 种 组 织 加 权 基 因 共 表达 网 络 ,结果 显示 埃及 伊 蚊 
这 9 种 组 织 涉及 11 个 基因 共 表 达 模 块 (图 2) ;将 转 
录 组 数据 对 应 的 不 同 组 织 ( 肉 尽 : 脑 、 触 角 、 咏 、 下 显 
须 、 卵 梨 ; 雄 蚊 : 脑 ` 触 角 、 前 足 . 中 足 \ 后 足 、 腹 部 末 
端 ) 作 为 研究 的 性 状 , 成 功 鉴 定 到 与 雌 蚊 触角 、 吃 、 
下 显 须 .卵巢 以 及 雄 蚊 脑 、 腹 部 末端 相关 的 6 个 组 织 
特异 性 模块 ,分 别 为 green, purple, red, turquoise, 
blue 和 magenta 模块 (图 2)。 蚊 虫 分 化 成 熟 的 多 个 
组 织 共 同 支 撑 其 追寻 寄主 吸取 血液 、 交 配 产 卵 、 传 
播 疾 病 等 生命 活动 过 程 。 研 究 表明 ,蚊虫 可 利用 多 
种 化 学 和 物理 线索 ,如 宿主 释放 的 气味 .C0,、 周 围 
环境 的 化 学 信号 以 及 雌 蚊 翅 产 生 的 特定 频率 等 ,并 
据 此 定位 宿主 的 位 置 .选择 配偶 和 找寻 产 卯 地 点 
(Hartberg, 1971; Cator et al., 2009; Matthews et 
al., 2016) 。 触 角 、 喉 和 下 显 须 在 此 过 程 中 起 着 极其 
重要 的 作用 ,因此 它们 被 认为 是 蚊虫 防治 控制 的 重 
要 靶 标 之 一 ( 张 唱 唱 , 2019 ) 。 其 中 触角 是 蚊虫 最 主 
要 的 嗅觉 融 官 ,在 成 蚊 触 角 上 有 许多 被 称 为 感受 需 
的 感觉 结构 ,这 是 进行 化 学 感知 的 物理 位 置 ( 张 静 
等 , 2019) ; 叭 是 蚊虫 在 进食 过 程 中 处 理 味 觉 代 码 的 
要 头 部 附属 物 ,其 作为 探 针 到 达 宿 主 皮 肤 下 的 血 
管 ( Maekawa et al., 2011) 。 下 上 显 须 是 蚊虫 主要 的 化 
学 感觉 器 官 ,可 感知 CO0， 和 挥发 性 气味 物质 等 
( Grant et al., 1995) 。 脑 是 蚊虫 神经 调节 的 重要 场 
所 ,也 是 目前 蚊虫 生物 防 探 研究 的 热点 之 一 ,其 内 会 
有 多 种 神经 递 质 ,如 乙酰 胆 碱 .CABA 以 及 谷 氨 酸 等 
( Matthews et al., 2016) 。 在 脑 的 高 级 处 理 中 心 可 整 
合 嗅觉 和 其 他 感觉 模式 ,最 终 将 这 些 信 号 转化 为 行 
为 ,其 可 能 是 蚊虫 防 控 的 一 个 新 靶 标 ( Riabinina et 
al., 2016) 。 而 卵 梨 腹部 末端 是 与 蚊虫 生殖 细胞 发 














































































































育 、 繁 殖 相 关 的 重要 场所 ,有 研究 认为 可 通过 将 肉 蚊 
转化 为 无 害 的 雄 蚊 的 策略 而 实现 蚊虫 防 控 ,提示 调 
节 蚊 虫 性 别 分 化 过 程 可 能 是 当前 防 控 蚊虫 的 一 个 重 
要 突破 口 (Hall et al., 2015)。 对 6 个 组 织 特 异性 模 
块 内 基因 的 注释 分 析 发 现 ,与 肉 成 蚊 触 角 特 异性 相 
KHY green 模块 具有 气味 结合 和 嗅觉 受 体 活性 等 功 
能 ;与 峻 成 蚊 嗓 特异 性 相关 的 purple 模块 具有 丝 氨 
酸 型 肽 链 内 切 酶 活性 和 丝氨酸 水 解 酶 活性 等 功能 ; 
与 雌 成 蚊 下 频 须 正 相 关 的 red. 模块 具有 丝氨酸 型 肽 
链 内 切 酶 活性 和 碳水 化 合 物 结合 等 功能 ;与 雄 成 蚊 
脑 特异 性 相关 的 blue 模块 在 生物 学 过 程 调节 、 信 和 号 
转 导 神经 系统 过 程 等 生物 学 过 程 中 发 挥 主要 作用 ; 
与 肉 成 蚊 卵 梨 特异 性 相关 的 turquoise 模块 在 DNA 
复制 细胞 分 裂 、 细 胞 周期 等 生物 学 过 程 中 发 挥 主要 
作用 ;与 雄 成 蚊 腹 部 末端 特异 性 相关 的 magenta 模 
Hut X INIM .草酸 代谢 \ 有 机 酸 代 谢 等 生物 学 过 程 
中 发 挥 主要 作用 ( 表 1 和 2)。 其 均 能 与 对 应 组 织 功 
能 相 联 系 , 说 明 组 织 特 异性 模块 的 识别 可 靠 ,可 用 于 
后 续 相 关 信息 的 挖掘 。 

核心 基因 在 基因 表达 网 络 中 具有 高 的 连接 度 ， 
其 往往 是 重要 的 作用 靶 点 和 研究 热点 ,在 生物 学 过 
程 中 起 着 关键 作用 。 但 通过 Cytoscape 构建 基因 互 
作 网 络 后 , 较 多 的 node 和 edge 常会 掩盖 真正 的 核 
心 基 因 ,影响 筛选 过 程 。 因 此 本 研究 采用 Cytoscape 
软件 中 的 CytoHubba 插件 ,进一步 缩小 了 6 个 组 织 
特异 性 模块 中 核心 基因 的 搜索 范围 , 共 筛 选 出 47 个 
在 模块 内 可 能 起 重要 作用 的 hub 基因 ( 表 3) 。 这 些 
基因 大 多 可 在 黑 腹 果 蝇 Drosophila melanogaster 中 找 
到 对 应 直系 同 源 物 。green 模块 内 AAELOO5810, 
AAEL007669 , AAELO10214 和 AAELO13197 在 黑 腹 果 
晶 中 的 直系 同 源 基因 分 别 为 CC14085 Anidh, 
CCDCI151 和 CG4362 ;blue 模块 内 hub 基因 在 黑 腹 果 
晶 中 均 有 直系 同 源 物 ，turquoise 模块 内 
AAELO008595 , AAELOO0961 和 AAELOO3226 在 黑 腹 果 
晶 中 有 直系 同 源 基 因 , 4 Mael, Clos 和 
Cpr49Aa-Cpr65Az; red 模 块 内  AAELO07033, 
AAELO00311 和 AAELO13417 在 黑 腹 果 蝇 中 直系 同 源 
基因 分 别 为 P5cr-2, Fon 和 CG30280 ; magenta 模块 
内 AAELO06427, AAEL006429, AAEL017567 和 
AAELOO04944 在 黑 腹 果 晶 中 的 直系 同 源 基因 分 别 为 
CGI7571-thetaTry, | thetaTry-CGY7571, | Nep5 和 
Muc12Ea ( Matthews et al., 2016) , 47 个 hub 基因 
中 ,29 个 基因 为 已 知 功能 基因 。 如 与 雌 成 蚁 触角 特 
异性 相关 的 green 模块 中 , 共 筛 选 出 9 个 hub 基因 
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(36 3) ,其 中 AAELO10426 和 AAELOO1510 被 描述 为 
气味 受 体 家 族 基 因 (He et al., 2015; Matthews et al., 
2016); AAELO14089 属于 离子 受 体 家 族 基因 ; 
AAELO10214 编码 卷曲 螺旋 结构 域 蛋白 ,AAELOT2424 
编码 微 管 蛋 白 。 与 肉 成 蚊 唆 特 异性 相关 的 purple 模 
块 中 , de Yeh 6 ^8 hub 基因 ( 表 3), 其 中 
AAEL000011 和 AAEL000039 属于 离子 受 体 家 族 基因 
( Matthews et al., 2016) ; AAELO02600 , AAELO13033 
fil AAELO13032 被 描述 为 丝氨酸 蛋白 酶 家 族 基因 
( Bartholomay et al., 2004; Matthews et al., 2016), 
与 肉 成 蚊 下 里 须 特异 性 相关 的 red 模块 中 , deor e 
出 9 个 hub 基因 ( 表 3) ,其 中 44EL011624 BETZA 
酸 和 蛋白 酶 家 族 基 因 ,AAEL000044 属于 乌 氮 酸 脱 羧 酶 
家 族 基 因 ,4451013417 编码 血管 生成 素 相关 和 蛋白 
( Matthews et al., 2016) ; AAELO07033 编码 A- 凝 集 素 
Ait x MV. AE, AAELOO3405. 为 微 管 动态 蛋白 的 调节 因 
T ,AAEL006920 被 描述 为 G 蛋白 偶 联 受 体 家 族 基 
(Bartholomay et al., 2004; Matthews et al., 2016) 。 
筛选 到 的 hub 基因 功能 与 组 织 功能 相关 ,表明 筛选 
到 的 基因 较为 可 靠 ,但 有 的 hub 基因 的 报道 研究 其 
少 ,提示 后 续 可 进一步 通过 blast, RT-PCR 以 及 
RNA 干扰 等 手段 对 筛选 到 的 hub 基因 功能 进行 预 
测 和 检验 。Matthews 等 (2016) 认 为 2 个 与 氨 转 运 有 
关 的 基因 AAELOO7373 和 AAELOOS046 可 能 参与 埃 
及 伊 蚊 寻找 宿主 过 程 中 氨 的 检测 ,本 研究 WGCNA 
分 析 结 果 发 现 AAELOOT373 被 划分 到 green 模块 ,其 
在 上 肉 / 雄 蚊 触 角 中 高 表达 , AAEL008046 被 划分 到 
blue 模块 ,其 在 上 肉 / 雄 蚊 脑 中 高 表达 ,这 与 Matthews 
等 (2016) 报 道 相 符 , 此 外 ,本 研究 构建 了 以 这 两 个 
基因 为 核心 的 局 部 网 络 ,揭示 其 基因 网 络 关 系 (图 
5). AAELO000252. 被 报道 为 肉 蚊 富 集 表达 基 
( Matthews et al.，2016), 在 本 人 研究 中 被 划分 到 
turquoise 模块 , 其 在 肉 成 蚊 卵 集中 高 表达 ; 
AAEL010631 被 报道 为 雄 蚊 富 集 表达 基因 ( Matthews 
et al., 2016) ,被 划分 到 magenta 模块 ,其 在 雄 成 蚊 腹 
部 末端 中 高 表达 。 

本 研究 是 基于 WGCNA 方法 和 转录 组 数据 对 蚊 
虫 不 同 组 织 基因 共 表 达 模 式 的 首次 探索 。 构 建 了 埃 
Je pisc tse 9 种 组 织 加 权 基 因 共 表达 网 络 ,获得 了 
11 个 基因 共 表 达 模 块 , 其 中 ,6 个 模块 与 组 织 特 异性 
相关 。 完 成 了 组 织 特异 性 模块 内 基因 的 功能 富 集 分 
Vr ,发 现 这 些 模块 的 生物 学 功能 均 与 对 应 组 织 功 能 
密切 相关 。 此 外 ,鉴定 了 6 个 组 织 特异 性 模块 内 共 
47 个 hub 基因 ,揭示 了 其 基因 互 作 网 络 关 系 。 本 研 





















































究 结 果 可 为 蚊虫 基因 共 表 达 模 式 分 析 提 供 方法 基 
础 ,对 探究 蚊虫 不 同 组 织 特有 的 基因 资源 信息 以 及 
功能 基因 生物 信息 研究 有 参考 价值 。 
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